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The double perovskite material Sr2FeMo06 has the rare and desirable combination of a half- 
metallic ground state with 100% spin polarization and ferrimagnetic T c ~ 420K, well above room 
temperature. In this two-part paper, we present a comprehensive theoretical study of the magnetic 
and electronic properties of half metallic double perovskites. In this paper we present exact diago- 
nalization calculations of the "fast" Mo electronic degrees coupled to "slow" Fe core spin fluctuations 
treated by classical Monte Carlo techniques. From the temperature dependence of the spin-resolved 
density of states, we show that the electronic polarization at the chemical potential is proportional 
to magnetization as a function of temperature. We also consider the effects of disorder and show 
that excess Fe leaves the ground state half-metallic while anti-site disorder greatly reduces the po- 
larization. In a companion paper titled "Effective Spin Hamiltonian and Disorder Effects in Double 
Perovskite Ferrimagnets" , we derive an effective classical spin Hamiltonian that provides a new 
framework for understanding the magnetic properties of half-metallic double perovskites including 
the effects of disorder. Our results on the dependence of the spin polarization on temperature and 
disorder has important implications for spintronics. 



I. INTRODUCTION 

Materials with half-metallic ground states in which 
conduction electrons are fully spin polarized and have fer- 
romagnetic transition temperatures T c well above room 
temperature are very rare in nature. Only two families 
of materials, the double perovskites and the Heussler al- 
loys, have exhibited this special combination of proper- 
ties. Consequently, they hold the potential for tremen- 
dous advancements in the field of spintronics as spin in- 
jectors and tunneling magnetoresistance devices. 

Here, we focus on double perovskites which have gen- 
erated considerable interest due to their close connec- 
tions to ternary perovskites. The perovskite family is 
known to exhibit a wide variety of exotic properties in- 
cluding high-T c superconductivity, colossal magnetore- 
sistance and ferroelectricity. Moreover, since these mate- 
rials are derived from the same general family there is the 
potential to grow lattice-matched layered materials with 
different functional properties in each layer. Double per- 
ovskites with the general formula A2BB'06 is a compos- 
ite of two different ternary perovskites ABO3 and AB'03 
arranged in a 3D checker board pattern. The additional 
flexibility of choosing two different transition metal ions 
in double perovskites opens up many new avenues of ma- 
terial exploration, like the juxtaposition of strong spin- 
orbit coupling and strong interaction by combining 5d 
and 3d transition metals. Already the range of proper- 
ties span metals to band insulators, and multi-band Mott 
insulators, as well as ferromagnets, antifcrromagnets, fer- 
roelectrics, multiferroics, and spin liquids r^— 

The most widely studied double perovskite S^FeMoOe 
(SFMO) has a half metallic ground state with a ferromag- 
netic transition temperature T c w 420K (Ref. ,216:]) which 
is well above room temperature. In spite of having a 
complex chemical structure, SFMO is a simple system to 



understand from a theoretical point of view. In contrast 
to ferromagnets like iron, there is a clear separation of 
the localized (Fe spins) and itinerant degrees of freedom 
(originating from Mo). Unlike the manganites, SFMO 
has neither Jahn- Teller distortions nor competing mag- 
netic ground states. Finally, in contrast to dilute mag- 
netic semiconductors, disorder is not an essential aspect 
of the theoretical problem. Previous theoretical work on 
half-metallic double perovskites includes pioneering T=0 
electronic structure calculations^, model Hamiltonians 
analyzed using various mean-field theories^—, and two- 
dimensional (2D) simulations^. 

In this paper, hereafter referred to as paper I, and 
its companion paper titled "Effective Spin Hamiltonian 
and Disorder Effects in Double Perovskite Ferrimagnets" , 
hereafter referred to as paper II we expand on our 
recent worki^ on the magnetic and electronic properties 
of SFMO. Several new results as well as important details 
which were omitted in our Letter— are covered. Broadly, 
paper I focuses on the properties of the itinerant quan- 
tum electrons and their effect on magnetism, while paper 
II describes in detail the derivation of the effective spin 
Hamiltonian and the results obtained from it. 

The main results presented in paper I are: (1) Using a 
variational analysis, we obtain a phase diagram as a func- 
tion of the parameters in the Hamiltonian and show that 
for the specific parameters of SFMO, it is firmly in the 
ferrimagnetic phase. (2) We present, for the first time, 
the temperature dependence of the spin resolved den- 
sity of states. The electronic polarization decreases from 
100% at T=0 with increasing temperature and vanishes 
above T c . More importantly, it shows that the polariza- 
tion at the chemical potential is proportional to the core 
spin magnetic moment as a function of temperature. This 
result is crucial because it allows us to infer electronic 
properties from magnetic properties obtained from the 



effective spin Hamiltonian. (3) Finaiiy, we present the 
dependence of the polarization on disorder. In particu- 
lar, we show that for Fe rich systems, the ground state 
remains half-metallic while anti-site disorder rapidly re- 
duces polarization. 

Paper I is organized as follows. In Sec. [TT1 we de- 
scribe the generalized double exchange Hamiltonian used 
to study SFMO. In Sec. lIII| we present a variational anal- 
ysis that describes the dependence of ground state mag- 
netic properties on Hamiltonian parameters. For param- 
eters relevant to SFMO, the ground state is deep in the 
ferrimagnetic phase. The results of a perturbative spin 
wave analysis are discussed in Sec. IIVI We show that the 
spin stiffness, which sets the scale for magnetic T c , is 
two orders of magnitude smaller than the electronic en- 
ergy scale, allowing us to separate the itinerant degrees 
of freedom from the localized spins in the spirit of the 
Born-Oppenheimer approximation. The temperature de- 
pendence of the spin resolved density of states (DOS) is 
presented in Sec. |V] Finally, in Sec. I VII we present the 
effect of disorder on electronic polarization and conclude 
with some remarks about future directions in Sec. IVHI 



II. MODEL HAMILTONIAN 

SFMO can be well understood in terms of a generalized 
double exchange model£ ~ 10 ' 12 . The large Hund's coupling 
on Fe 3+ (3<i 5 ) leads to a local S — 5/2 core spin. Since 
locally all the spin up states on Fe are occupied, the only 
channel for the 4c? 1 electron on Mo 5+ to delocalize is by 
hopping from one Mo site to the next via the unoccupied 
Fe down states. This naturally leads to an antiferromag- 
netic coupling between the core Fe spins and the itinerant 
Mo electrons due to Pauli's exclusion principle. The con- 
duction band is formed by hybridization of the Fe t2 g \ r 
and Mo ti g orbitals via oxygen. Symmetry dictates that 
d a /3 electrons delocalize only in the (a, j3) planed where 
a, /3 = x, y, z. The model Hamiltonian describing SFMO 
is 

H = ~t ^ (eiad^Cja + h.c.) 

E c W° +A T, d l d u (!) 

where di (cj) denotes the fermion operator on the i th Fe 
(Mo) site. The electronic spin on Fe site is quantized 
along the direction of the local spin, whereas on Mo site 
the quantization is along a global z-axis. t is the near- 
est neighbor Fe-Mo hopping amplitude, t' is the direct 
hopping amplitude between Mo sites and A is the charge 
transfer energy between Fe t2 g \. and Mo ti g states. In the 
global frame of the Mo spins, the orientation of the i th Fe 
core spin is given by (ftj, <pi) and it determines the effec- 
tive hopping amplitude between Fe and Mo sites through 
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Fe(3d 5 ) Mo(4d 1 ) Fe Mo 

FIG. 1: Schematic showing energy levels at the transition 
metal sites in two unit cells (formula units) of SFMO. The Fe 
sites have localized 5 = 5/2 core spins (red arrows), treated 
as classical vectors with orientation (9,<j>). The up and down 
sectors on the Fe site are split by a combination of the Hund's 
coupling Jh and onsite Hubbard term U . The parameters t, 
t' , and A of the Hamiltonian ([T]l. governing the dynamics of 
the itinerant electrons (blue arrows) in t^g orbitals, are also 
shown. 

eif and Cii which are defined as 

e*f = -sin(0i/2)exp(i<&/2) 

e 4 = cos(0i/2)exp(-*&/2) (2) 

A schematic of the level structure is shown in Fig. [TJ In 
pure SFMO, we ignore direct Fe-Fe hopping and Fe-Fe 
superexchange because the Fe sites are far apartii and 
the spatial extent of the 3d orbitals is much smaller com- 
pared to 4d orbitals. However, these can be important 
in the presence of disorder since two Fe sites can be right 
next to each other. 

At first glance, the double exchange model for double 
perovskites looks like the antiferromagnetic Kondo lat- 
tice model 15 . However, there are significant differences: 
In the case of the antiferromagnetic Kondo lattice, the 
local moments are quantum degrees of freedom whereas 
in SFMO the local moments are large (S — 5/2) and 
treated classically. Another consequence of the large lo- 
cal moments in SFMO is that they cannot be completely 
screened by the available conduction electrons. Finally, 
as mentioned above, Pauli's exclusion principle is respon- 
sible for the antiferromagnetic coupling between the local 
moments and the itinerant electrons in SFMO. In con- 
trast, the coupling in Kondo lattice arises from antifer- 
romagnetic exchange interaction and is typically small 
compared to the band width. 

III. VARIATIONAL ANALYSIS 

We begin by studying the T=0 properties of the quan- 
tum Hamiltonian in Eq. ([T]). As a function of Hamilto- 
nian parameters, we explore the relative stability of var- 
ious magnetic phases shown in Fig. [21c): ferromagnetic 
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FIG. 2: Variational calculation: (a) A-n magnetic phase di- 
agram for t'=0; here A is the charge transfer offset between 
the Mo tig orbitals and the Fe t2 g i orbitals, and n is the fill- 
ing, (b) t'-n phase diagram for A=0. Solid lines are first 
order phase boundaries and the black dot indicates SFMO 
parameters. Schematic of the variational magnetic phases, 
ferromagnet (FM), stripe antiferromagnet (AF1), and usual 
Neel antiferromagnet (AF2) are shown in (c). The red arrows 
indicate Fe spins and blue dots indicate Mo sites. 



(FM), stripe antiferromagnet (AF1) and Neel antiferrmo- 
magnet (AF2). Note that the FM phase has induced 
moments on Mo sites which are aligned in the opposite 
direction. So, technically, it is a ferrimagnetic state, but 
we will focus on the core spins and refer to this state 
as FM. We use the same nomenclature in paper II also. 
In Fig. |U(a), we present the A — n phase diagram for 
t' = 0, where < n < 3 is the electron filling. SFMO 
corresponds to n = 1/3. Part (b) shows the t'-n phase 
diagram for A=0. The phase boundaries are determined 
by the relative energy of conduction electrons in the var- 
ious variational spin backgrounds. For the FM state, 
the energy is calculated by integrating over the filled 
states in the band structure shown in Fig. U{a). Details 
for the AF1 and AF2 phases are provided in Appendix 
[Al Our result in Fig. Ufa) is consistent with previous 
calculations^. 

The scale for the magnetic T c is set by t, and we will 
show in paper II that choosing t — 0.27 eV, consistent 
with band structure calculations 6 -, leads to the experi- 
mental T c « 420K of SFMO. For the other Hamiltonian 
parameters, we use i'/i = 0.1 and A/t = 2.5 which are 
also in agreement with band structure calculations 6 -. It is 
clearly seen from Fig. [U[a) and (b) that, for the param- 
eters relevant to SFMO (indicated by black dots), the 
ground state is deep inside the FM phase. This justifies 
our claim that SFMO does not have competing magnetic 
phases. 

It is, however, interesting to note that t' has a ma- 
jor role in determining the phase boundaries. This is 
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FIG. 3: Comparison of spin wave dispersion from perturba- 
tion theory [Eq. Q] with that obtained from exact diagonal- 
ization of Eq. JTJ) in a spin wave background. We find excellent 
agreement with no fitting parameters. Also note that the en- 
ergy scale of magnetic interactions is two orders of magnitude 
smaller than the electronic band width (w 8t). 



due to the fact that the effect of t' is very different in 
the different phases. For t' = 0, the AF1 phase has 
one dimensional bands and increasing \t'\ introduces two 
dimensional hopping which changes the nature of the 
bands dramatically. However, in the case of FM and AF2 
phases, the bands are two dimensional in nature even for 
if = and, therefore, the effect of t' is not as strong as in 
the AF1 phase. The t' < region of our phase diagram 
can be mapped on to the t' > region by the symmetry 
of the Hamiltonian: E(t' , A, n) = E(-t' , -A, 3 - n). 

On the other hand, A has very little effect on the T=0 
phase diagram because its effect is very similar in all the 
phases. We, however, expect increasing A to significantly 
reduce T c . These phase diagrams are useful in guiding 
materials search with optimized parameters. 



IV. SPIN WAVES 

In this section, we calculate the spin wave dispersion 
and the spin stiffness of SFMO using a perturbative anal- 
ysis. At low temperatures, the core Fe spins fluctuate 
about the fully magnetized ferromagnetic state and these 
fluctuations affect the mobile electrons. To lowest order, 
we can view the core spin fluctuations as generating spin 
wave configurations that are static on the time scale of 
the electronic degrees of freedom. This separation of time 
scales for the core spins and the itinerant electrons will 
be justified a posteriori, and plays an important role for 
later results. 

The classical core spins in a frozen spin wave can be 
described as 
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(3) 



where q is the wave vector of the spin wave. We assume 
that the angle 9 with respect to the quantization axis of 
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FIG. 4: (a) Electronic polarized band structure of FM ground state (blue indicates spin down and red, spin up). At SFMO 
filling, only the lowest spin-down band is occupied, thereby confirming the half-metallic ground state, (b), (c) and (d) shows the 
spin resolved density of states at T=0, T~T c /2 and T>T c respectively. The black horizontal line indicates chemical potential 
at SFMO filling. 



the FM ground state is small and we explicitly calculate 
the corrections to the energy up to 0(9 2 ). In Eq. ([T]) the 
fermion operators on the Fe sites are described with re- 
spect to the local quantization axis. In the analysis here, 
for small angular perturbations of the Fe spins, it is con- 
venient to redefine the creation(annihilation) operators 
on Fe sites in the same global frame as the operators on 
Mo sites. We choose 

dl = cos (9/2) fl - sin {6/2) (4) 

where ^(f^) is the creation(annihilation) operator in 
the global frame for an electron with spin a on the i th Fe 
ion. The Hamiltonian in Eq. ([TJ can now be rewritten in 
terms of these new operators. Keeping only terms up to 
0(0 2 ), we get 

H = H a + 0H 1 + 6 2 H 2 (5) 

where 

E • h - c ) ( 6 ) 

H ^ = \ E ( e ~ iq ' r, 4r^ + ^'%^t + h.c) (7) 

(i,3) 

^ = lE( f 4 C ^- f ^T + H (8) 

Here Hq is the unperturbed Hamiltonian, while Hi 
describes the hybridization of spin up and spin down 
orbitals which is unique to the double perovskites. Fi- 
nally, Hi contains terms responsible for narrowing the 
spin down conduction band while allowing spin up elec- 
trons to delocalize. The details of the calculation are 
given in Appendix [B] 



One subtlety that is worth pointing out is that we are 
now working in an over complete basis. The Pauli block- 
ing of local spin up states on Fe sites is enforced in the 
global frame by constraining fL and f,^ operators to ap- 
pear only in specific linear combinations that correspond 
to the local down spin operators. 

We obtain the spin wave dispersion by taking the sec- 
ond derivative with respect to 9 of the total change in 
energy: 

k|</c F 

where the sum is taken over all occupied levels. <5e(k, q) 
[see Eq. (|B14I) ] is the change in energy up to 0(9 2 ) of the 
state in the conduction band labeled by k. The resulting 
spin wave dispersion is shown in Fig. [3] It agrees very 
well with the spin wave dispersion calculated by exact 
diagonalization on a finite lattice, also shown in Fig. [3] 
No fitting parameters are used. The analytical result 
from the perturbative analysis, presented here, has the 
advantage of being not restricted by finite system sizes. 
In paper II, we use the spin wave dispersion to extract 
the parameters of the effective spin Hamiltonian. 

We have further calculated the spin stiffness, J g = 
lim q ^o(3 2 -E(q)/9g 2 ). For parameters relevant to SFMO, 
we find that J e g = — 0.035i which is about two orders of 
magnitude smaller than the band width (W ~ 8t) of the 
itinerant electrons. This justifies our initial assumption 
that on the time scale of the electrons, the spin waves 
can be approximated by static spin configurations. 

V. TEMPERATURE DEPENDENT DENSITY 
OF STATES 

We use a method that combines exact diagonalization 
with Monte Carlo (ED-f-MC) to calculate temperature 




FIG. 5: (a) Core spin magnetization (M) and polarization 
of the conduction electrons at the chemical potential (P) of 
SFMO as a function of temperature calculated using ED+MC 
method on an 8 x 8 lattice, (b) Parametric plot of P (T) against 
normalized M(T). It is clear that P(T) is proportional to 
M(T). The dashed line indicates exact proportionality. 

dependent properties of SFMO 1 -^. For each spin con- 
figuration, the electronic energy is calculated by exact 
diagonalization which is then used to update the spin 
configuration using Monte Carlo algorithm. The assump- 
tion that the fast electrons relax immediately to the given 
spin texture has already been justified in Sec. IIVI by the 
clear separation of time scales for the local and itinerant 
degrees of freedom. At each Monte Carlo step, a new 
random spin orientation is generated using Marsaglia's 
methoo— i and acceptance is based on the metropolis algo- 
rithm. All calculations are done on lattices up to 16x16, 
and twisted boundary conditions are used to minimize 
finite-size effects. 

We use ED+MC method to calculate, for the first time, 
the spin resolved density of states (DOS) as a function of 
temperature, shown in Fig. [4] At all temperatures, the 
spin quantization axis is defined along the direction of 
magnetization, which is the natural axis for a ferromag- 
net. 

Fig. QJa) shows the spin polarized bands of the FM 
ground state (red indicates spin down while blue denotes 
spin up). The bonding band of the Fe ti g \, and Mo t2 g \ r 
orbitals form the conduction band while the anti-bonding 
band is pushed up in energy. The spin up band in the 
middle comes from the Mo t2 g \ orbitals and they don't 
hybridize with the spin down orbitals in the perfect FM 
state. 

We find that for T=0, only the spin down bonding 
band is occupied and SFMO is a half metal, in agreement 
with photoemission experiment^ and electronic struc- 
ture calculations-? 6 . For 0<T<T c , the broken time re- 
versal symmetry leads to very different DOS for spin up 
and spin down sectors with the DOS at chemical poten- 
tial dominated by spin down. In clear distinction from 
the strictly T=0 case, both spin sectors have non-zero 
DOS at all energies. Finally, for T>T C , there is no pre- 
ferred spin direction and DOS for spin up and down are 
identical. 

As seen in Fig. |U(b),(c) and (d), the DOS varies 




FIG. 6: (a) Zero temperature polarization of conduction elec- 
trons, P(0), in off-stoichiometric SFMO with general formula 
Laa;Sr2-xFei+ a Moi- a 06 (y > is Fe rich while y < is Mo 
rich), (b) P(0) as a function of anti-site disorder. Anti-site 
disorder is characterized by S which is the fraction of Fe on 
the wrong sublattice. Fe rich systems remain half-metallic 
while anti-site disorder or Mo excess rapidly decreases the 
polarization. 

smoothly with temperature. This has an important con- 
sequence that the polarization of the conduction elec- 
trons at the chemical potential, P=(N^-N-j-)/(N^+N-j-) 
where N CT is the density of states of spin a at the chemical 
potential, is proportional to the magnetization of Fe core 
spins, M, as a function of temperature^. In Fig.[SJa), M 
and P arc plotted as a function of temperature. For bet- 
ter visualization of the proportionality, we have shown a 
parametric plot of P against a normalized M in Fig.[5](b), 
with T as the implicit parameter. A linear fit describes 
the data very well. The proportionality of P and M 
is crucial experimentally because P(T) is the quantity 
of interest in spintronics applications but is difficult to 
measure. Our result allows direct inference of polariza- 
tion from the magnetic moment, a much simpler quantity 
to measure experimentally. From a theoretical point of 
view also, the proportionality of P and M allows us to 
focus only on the magnetism. In paper II, we derive an 
effective classical spin Hamiltonian which describes the 
thermodynamics of the Fe core spins. It facilitates accu- 
rate calculation of magnetic properties and, by virtue of 
the proportionality, also provides realistic results for the 
electronic polarization. 



VI. EFFECT OF DISORDER 

In this section, we discuss the effect of disorder on elec- 
tronic polarization at low T. Its effect on magnetization 
is described in paper II. We consider three types of disor- 
der: 1) excess Fe, 2) excess Mo, and 3) anti-site disorder. 

Excess Fe: The general formula of off-stoichiometric 
SFMO is La a; Sr2_ 2 ;Fei+j ) Moi_j,06 and Fe rich systems 
correspond to y > 0. Some Mo is substituted by Fe, 
which generates antiferromagnetic superexchange cou- 
pling S(S + 1)Jse ~ 34meV between Fe on neighbor- 
ing sites. The exchange coupling J$e is estimated from 
that of the AF insulator LaFeOs with Tn ~ 750K using 
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S(S + 1)Jse = Since Mo is also the source 

of itinerant electrons, excess Fe decreases the filling. In 
Fig.|UJa) we show the electronic polarization at zero tem- 
perature, P(0), for Fe rich SFMO. An important ob- 
servation is that the conduction electrons remains fully 
spin polarized. The persistent half-metallicity can be un- 
derstood intuitively from the fact that at low tempera- 
tures the strong superexchange interaction locks the ex- 
cess Fe on the wrong sublattice in a perfect antiferro- 
magnetic configuration with respect to its neighboring 
Fe sites. Consequently, the extra Fe sites are prevented 
from participating in the derealization of itinerant elec- 
trons, while the rest of the lattice continues to have a 
ferromagnetic ground state with spin polarized conduc- 
tion electrons. 

Excess Mo: When Mo substitutes for Fe, there are re- 
gions with Fe spins that are much farther apart than in 
the perfect lattice. In addition, the density of carriers 
in Mo-rich systems is higher than that of pure SFMO. 
Fig. [Ma) shows that unlike excess Fe, P(0) rapidly de- 
creases as a function of excess Mo (y < 0), consistent 
with DFT calculations^. The Mo rich regions constitute 
small regions that can be described by tight binding lat- 
tices with no preferred spin direction. They introduce 
states with both spins in the entire energy range of the 
conduction band which therefore reduces the polarization 
at the chemical potential. 

Anti-site: Finally, anti-site disorder, which is the most 
common type of disorder, arises when some Fe sites ex- 
change positions with Mo sites. It can be thought of as 
a combination of Fe rich and Mo rich regions while keep- 
ing the overall stoichiometry unchanged. Since anti-site 
disorder introduces Mo rich regions, P(0) decreases with 
increasing anti-site disorder as shown in Fig. |5fb). The 
extent of anti-site disorder is parametrized by S which 
is defined as the fraction of Fe on the wrong sublattice; 
complete disorder corresponds to S = 0.5. While anti-site 
disorder behaves quantitatively like excess Mo, the elec- 
tronic polarization decreases much faster with anti-site 
disorder. Given that fact that it is also the most common 
form of disorder, it is of crucial importance that anti-site 
disorder be minimized in order to get the high electronic 
polarization required for spintronics applications. 



VII. CONCLUSION 

In conclusion, we have presented here a general frame- 
work for understanding half-metallic double perovskites. 
While the generalized double exchange model predicts 
other magnetic phases as a function of Hamiltonian pa- 
rameters, we have shown that for SFMO only the FM 
phase is relevant. We have calculated for the first time 
the temperature dependence of the spin resolved DOS 
and found a proportionality between the temperature de- 
pendence of the electronic polarization and the magneti- 
zation which is a significant result. It offers a much sim- 
pler method for determining the polarization. Finally, 



we have shown that Fe rich systems have a half-metallic 
ground state while anti-site disorder greatly reduces the 
polarization. Such understanding is crucial for spintron- 
ics applications. 

The results of paper I become the starting point for 
paper II. Motivated by the proportionality between the 
electronic polarization and the core spin magnetic mo- 
ment, in paper II, we focus entirely on the large local 
spins on Fe sites and infer electronic properties from 
the magnetization. We derive an the effective magnetic 
Hamiltonian describing the thermodynamics of the clas- 
sical spins. The effective Hamiltonian offers a new frame- 
work for understanding the magnetic properties in half 
metallic double perovskites, including the effects of dis- 
order on the saturation value of magnetization and the 
T c . We also take advantage of the fact that Fe excess 
does not change the half metallic ground state to pro- 
pose a novel way of increasing T c without sacrificing the 
polarization. 
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Appendix A: Variational Analysis 

The electronic bands for the variational magnetic 
states can be obtained by Fourier transforming the dou- 
ble exchange Hamiltonian in Eq. (|TJ) into momentum 
space. The minimum dimension of the Hamiltonian ma- 
trix in momentum space is related with the periodicity 
of the configuration. FM has periodicity of one unit cell 
and every unit cell has three states, therefore Hp Af (k) is 
a 3x3 matrix. However AFI and AF2 has a minimum 
periodicity of two unit cell. Thus HafiO*) and Haf2{^) 
are six dimensional. The Hamiltonian matrices are 

/ A -2tsi(k)' 

H FM (k) = -2t'g 2 (k) 

\-2*<?i(k) -2t'g 2 (k) / 



where 



3i(k) = cos I — — I +COS ' 



32 (k) = cos(k x a) + cos(k y a) 
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FIG. 7: J c ff as a function of Hamiltonian parameters: (a) J e fr 
vs t' for A = 2M, (b) J cff vs A for t' = O.lt. Filling is fixed 
at SFMO value (n = 1/3). 



= -2tcos((k x + k y )a/2) 
Vk = —2tcos((k x — k y )a/2) 
Afe = — 2t (cos(k x a) + cos(fcj,a)) 



(A3) 



The relevant energy bands are obtained by diagonaliz- 
ing these matrices as a function of k, and the energy of 
the conduction electrons is calculated by integrating over 
filled levels. The relative energy of FM, AF1 and AF2 
determines the phase boundaries in Fig. fj] 



Appendix B: Second order perturbation theory 

We describe here the perturbative calculation of the 
spin wave dispersion and the spin stiffness of SFMO. The 
Fe core spins cant by a small angle 9 in a spin wave 
configuration [Eq. (|3j|], and we calculate corrections to 
the energy levels of the FM ground state up to 0(9 2 ). 
The unperturbed Hamiltonian Hq in Eq. ([6]) gives four 
bands with eigenvectors given by 



ot(k) 

4(k) 
4(k) 
4(k) 



a(k)f£j, + 
-/3(k)f^ 



+■ «( k ) c U 



f f 



(Bl) 



The corresponding eigenvalues are 



ca(k) 
e 4 (k) 



A 
~2 

^-t'. 92 (k) + r(k) 



i' 52 (k) - r(k) 
t'saO 

2i' 52 (k) 



where 
a(k) 

(3(k) 

r(k) 

A(k) 
5i (k) 
92 (k) 



e 3 (k)-ei(k) 
2tgi(k) 

V ^k) 

V /4^g2(k) + (A/2 + t' ff2 (k)) 2 
( ei (k)-e 3 (k)) 2 + (2i 5l (k)) 2 



cos 



cos 



kyd 



cos(k x a) + cos(k y a) 



(B3) 

(B4) 

(B5) 
(B6) 

(B7) 

(B8) 



The first band is the bonding Fe^-Moj, band. The second 
is the anti-bonding band. The third and fourth are the 
Mo| and Fe^ bands respectively. The distance between 
two Fe ions is a. For SFMO filling only the lowest band 
ei(k) is occupied. 

We next describe how the energy of the lowest band 
ei(k) is affected by H\ and 7?2- First order correction in 
the canting angle 9 is 



e^k.q) = 0<ai(k)|ffi|ai(k)> =0 



(B9) 



The second order correction has several contributions. 
One of them is 



<ai(k)|ff 2 |ai(k)> =i 5 i(k)a(k)/3(k) 



(B10) 



The mixing of the lowest band with the anti-bonding 
band gives 



v l( Q2 (kQ|g 1 | Ql (k))| 2 
% ei(k)-c a (kO 

The mixing with the MoJ. band gives 

y, |<a 3 (k')|ffi|ai(k))| 2 _ t 2 . 9 2 (k-q)a 2 (k) 

k' 



A 



(B2) 



ei (k)-e 3 (k') ei(k)-ea(k-q) 

Finally, the mixing with Fej" band gives 

v |(a 4 (kO|g 1 [a 1 (k))| 2 = t 2 5 2 (k)/3 2 (k) 
^ ex(k)-e 4 (kO 



ei(k)-A 



(Bll) 



(B12) 



(B13) 
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Upon collecting all the second order correction terms and 
simplifying them algebraically, we get the energy correc- 
tion in lowest band up to O(0 2 ) 

5e(k,q)=e( 1) (k,q) + ef ) (k,q) 



A(k) 



(e 3 (k)- ei (k)) 



gi 2 (k-q) ( £l (k)-6 3 (k)) 2 
A(k) ei (k)-e 3 (k-q) 



(B14) 



Notice that the energy correction in Eq. (|B14|) has two 
terms. The first term comes from narrowing of conduc- 



tion band in the spin wave background. It increases spin 
stiffness. The second term comes from hybridization of 
spin down conduction band with the Mo-j- band and it 
reduces the spin stiffness. As shown in Eq. ©, the spin 
wave dispersion can be calculated by summing Se(k, q) 
over all filled states. 

We can also use the result from perturbative anal- 
ysis to calculate change in spin stiffness, J c g — 
lim q _ > o(5 2 i?(q)/9(7 2 ), as a function of Hamiltonian pa- 
rameters. Fig. [TJa) shows the dependence of J e g on t' 
while Fig. EJb) shows how J e g changes with A. Increas- 
ing either t 1 or A decreases J c g as expected. 
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